Extension of the Fluctuation-Dissipation theorem to the physical aging of a model 
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We present evidence in favor of the possibility of treating 
an out-of-equilibrium supercooled simple liquid as a system in 
quasi-equilibrium. Two different temperatures, one controlled 
by the external bath and one internally selected by the sys- 
tem characterize the quasi-equilibrium state. The value of 
the internal temperature is explicitly calculated within the 
inherent structure thermodynamic formalism. We find that 
the internal temperature controls the relation between the re- 
sponse to an external perturbation and the long-time decay 
of fluctuations in the liquid. 

In the last decade, several efforts have been devoted 
to the understanding of the glass-transition phenomenon 
, one of the open fundamental problems of con- 
densed matter. In particular, the thermodynamics of 
supercooled liquids and the relations between dynamics 
and thermodynamics have received considerable interest. 
Novel approaches and detailed analyses of computer 
generated configurations of supercooled liquids |£-13 are 
supporting and formalizing the picture — rooted in old 
ideas jL4| — of a glass as a system trapped in one lo- 
cal free energy basin. Correspondingly, the equilibrium 
dynamics of deep supercooled liquids is interpreted as 
motion among different basins. 

In this Letter, we show that the thermodynamic ap- 
proach also allows us to interpret the behavior of super- 
cooled model liquids under non-equilibrium (aging) situ- 
ations. The model liquid we study is the Lennard Jones 
binary (80:20) mixture (BMLJ) which has been exten- 
sively investigated in the past The formalism 
we adopt is the so-called inherent structure (IS) formal- 
ism J!?]] , which focus on the local minima of the potential 
energy surface and on the corresponding basins of attrac- 
tion. In this formalism, each point in the 3A-dimensional 
configuration space is unambiguously associated to the 
basin of attraction of the closest IS. The configuration 
space is partitioned in a sum of basins (usually labeled 
according to the value ejs of the potential energy in the 
IS) and the partition function is written as a sum of 
basin partition functions. In the thermodynamic limit, 
the system free energy F(T) is written as |17| , [l8]] 

F(T) = -TS conf (e IS (T)) + f basin (T, e IS (T)) (1) 

where —TS con f(eis(T)) account for the entropic contri- 
bution arising from the number of basins of depth e/s 
and fbasin(T,eis) describes the free energy of the sys- 
tem constrained in one characteristic e/s basin. The 
separation of the free energy in two parts reflects the 



separation of time scales which is observed in the su- 
percooled liquid state. Indeed, the intra-basin dynamics 
is much faster (ps-time scale in real liquids) compared 
to the inter-basin dynamics (whose time scale diverges 
upon cooling). We stress that the IS expression for the 
liquid free energy, Eq. (0) is analogous to the free-energy 
derived using the Thouless-Anderson-Palmer (TAP) JlSf] 
approximation in the p— spin models, once a basin is iden- 
tified with one TAP solution and fbasin(T, e/s) with the 
TAP free-energy. 

In the case of glass-forming liquids, the intrabasin free- 
energy f basin is usually written as e/s + fvtb(T,e IS ) to 
separate the value of the potential energy in the minimum 
(e/s) from the vibrational free-energy contribution. At 
low T, fvib(T,ejs) can be calculated in the harmonic 
approximation, expanding the potential energy around 
the IS configuration, as 



3JV-3 



fvib(T,eis) = k B T ln(hu>i(ei S )/k B T) 

i=i 



(2) 



where u>i is the frequency of the z-th normal mode p(| . If 
all basins had the same curvature, then f v a would be only 
function of T. In the BMLJ system, basins with different 
depth have different curvatures and hence f v n, is a func- 
tion of both T and e/s ■ Fig. |] shows the T-dependence 
of fbasin(T, eis) in the harmonic approximation, for dif- 
ferent e/s values. Curves for different e/s values are not 
parallel, since the basin curvatures depend on the basin 
depth, as shown in the inset of the same figure. 

In equilibrium, at each temperature T eq , the system 
populates basins of different depth eis(T eq ) The 
eis(T eq ) value is fixed by the condition of F being a min- 
imum, i.e. by 



dF 



-T, 



dS con f(eis) dfbasin(T eq ,eis) 
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(3) 



where the term ( 9 g b e ° 3 s '" ) describes the free energy change 
associated to the change in the basin shape with e/s and 



the term ( 



us. 



de IS 



describes the entropic change associ- 



ated to the different basin degeneracy. As shown in Ref. 

the onset of slow dynamics correlates with a sharp 
drop in the T-dependence of e/s- This observation holds 
for all model systems studied so far. For the BMLJ sys- 
tem, the e/s basins which are populated at the different 
temperatures (in the T-range accessible to equilibrium 
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simulations) are indicated in Fig. [|. As indicated by 
Eq. (||), the T-dependence of eis arises from a balance 
between the change of basin free energy and the change 
in configurational entropy. Of course, the monotonic re- 
lation eis(T eq ) can be inverted to give the temperature at 
which the equilibrium system populates basins of depth 
eis, T eq (e IS ) f22|. 

The analysis of the non-equilibrium dynamics pre- 
sented in this Letter is based on the assumption that 
the separation of intrabasin and interbasin time scales — 
which characterize equilibrium supercooled liquid states 
- retains its validity during aging. In particular, we 
assume that after a temperature jump from Ti to Tf, 
the vibrational intrabasin dynamics thermalizes quickly 
to the thermostat value Tf. The thermalization of the 
entire system is instead very slow, requiring a search for 
the low energy minima. Such a search in configuration 
space, which is the essence of the aging phenomenon in 
liquids, is so slow that quasi-equilibration on basins of 
depth e/5 might be faster than the decrease of the eis 
value. 

In the case of BMLJ, the out-of-equilibrium dy- 
namics following a T-jump have been recently studied 
p3| , p2]| . Ref . indeed suggested the possibility that the 
equilibration in configuration space proceeds via quasi- 
equilibrium steps. The proposed equilibration process is 
schematically depitched with arrows in Fig. 51 for the 
case of a jump from Tj = 0.8 to Tf = 0.25. In a time 
much shorter than the any basin change, fbasin assumes 
the value characteristic of the final temperature (full line 
arrow in Fig. [l]). The fast equilibration of the intra-basin 
degrees of freedom is then followed by a much slower pro- 
cess (dashed arrow in Fig. |l|) during which the system 
populates deeper and deeper eis levels. For the BMLJ 
case, eis(t) after a T-jump is reported in Fig. [|. 

If the hypothesis of quasi-equilibrium is correct, then 
we can ask which is the value of the internal temperature 
Ti n t(eis, Tf) selected by the system when it is populating 
basins of depth ejg 0. To calculate T int (eis,Tf), we 
again search for solutions of Eq. (^) but, in contrast to 
the equilibrium case, we consider the value eis to be fixed 
and solve for the unknown temperature, obtaining 



Ti n t(eis, Tf) 



1 + g§77fvib(Tf,e IS ) 



<:)<:■: j 



;S con f(eis) 



(4) 



Note that, differently from Eq. ([|), is now evalu- 
ated at the thermostat temperature Tf, since the fast 
intrabasin degrees of freedom are already thermalized to 
Tf. This expression for Ti„ t coincides with the expres- 
sion proposed by Franz and Virasoro p5[ in the context 
of p— spin systems, once the basin free energy is identified 
with the TAP free energy. We note that ^° nf can be 
evaluated from equilibrium conditions Eq. (ffl, and thus, 
once a model for is chosen, Ti nt (eis,Tf) can be cal- 
culated. If Tf is small (as it usually is), the harmonic 



approximation for f v n, can be confidently used. In this 
case, from Eq. fl) and Eq. (H) we obtain 



T mt {eis,Tf) 



i+k B T f j:, 



dln[hui/k B Tf] 



1 + k B T eq (ei S ) J2: 



dln[hi 



de IS 



(5) 



Fig. g-left shows T int (e IS , Tf) for the BMLJ system. We 
note that, if basin curvatures were independent on the 
eis value, then the derivatives in Eq. (^|) would be zero. 
Hence, Ti nt (eis,Tf) would not depend of Tf and Ti nt 
would coincide with T eq (eis). In this limit, Tint plays 
the same role as the Active temperature introduced in 
the analysis of experimental data in aging systems [^6j . 
We note in passing that an important by-product of the 
present approach is a free-energy expression for an out- 
of-equilibrium liquid, that depends only on eis and Tf. 
Such expression offers a detailed interpretation of the free 
energy expression for glassy systems proposed in recent 
years R], 

To test the predictions of Eq. (||j|) — i.e. the assump- 
tion of quasi-equilibrium in the aging system — we study 
the response of the liquid, described by an Hamiltonian 
Hq, to an external perturbation switched on at t = t w . If 
the perturbation adds the term Hp — —VoB{r N )0{t — t w ) 
to the Hamiltonian (where 0(t) is the Heaviside step func- 
tion), linear response theory predicts that the time evo- 
lution of any variable A(r N ) conjugated to B is given by 



(A{r)) 



Vo 

k R T 



[<A(r)B(0)) - <A(0)B(0)>]„] (6) 



where r = t — t w , (• • •) is the ensemble average over the 
perturbed system (H + Hp) and (• • -) is the ensem- 
ble average over the unperturbed system (Hq). We chose 
B = {Pi + PT), where p£ = £f Q e^< / yf{N) is the 
Fourier transform component of the density of a parti- 
cles at wavevector k, and study the response of A = p£ . 
In this case, (A(t)B(0))o coincides with the dynamical 
structure factor Sg a (t) = (Pk WPk*(°))o- Thus, with the 
present choice of A and B, linear response theory predicts 



<p£(r)> = -t%5£°('-) " SnO)] (7) 

Eq. (0) — also referred as fluctuation-dissipation theo- 
rem — is particularly relevant for our purposes, since 
it predicts that the response of the system is propor- 
tional to T~ l and thus offers an independent way to con- 
firm the validity of the quasi-equilibrium hypothesis. In- 
deed, the quasi-equilibrium hypothesis predicts that for 
short times (i.e. in the time region where the correla- 
tion function assumes values between the value in zero 
and the plateau value) the relation between correlation 
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and response, Eq. (Q), should be controlled by Vo/ksTf, 
since the intrabasin dynamics is probed. Similarly, for 
long times (i.e. in the time region where the correla- 
tion function assumes values smaller than the plateau 
value) the relation should be controlled by V /kBT int , 
since the inter-basin dynamics is now probed. Thus, we 
predict that switching on the perturbing field when the 
system is populating basins of depth e/s, the calculated 
Ti n t should coincide with the temperature at which the 
system responds to the external perturbation for times 
longer than the vibrational dynamics. 

Fig. ||-left shows Sg a (r) and p£(t) for two different t w , 
for the BMLJ case, with T % = 0.8 and T f = 0.25. The 
reported data are averaged over 300 different quench- 
realizations and over more than 60 different independent 
perturbations Hp (but with the same k modulus) for 
each quench-realization (28). Both S£ a (r) and p^(r) 
show the two-step relaxation characteristic of the super- 
cooled state, which has been associated to the separation 
of intra-basin and inter-basin motion. Fig. fright shows 
the corresponding response vs. correlation plots, Eq. (Q). 
At short time (intra-basin motion) pk(i) vs. Sk (t) is lin- 
ear with the expected TJ slope, properly describing the 
equilibrium condition of the vibrational dynamics with 
the external reservoir. At larger times, the intra basin 
motion sets-in and indeed the slope of p^(r) vs. S^ a (r) 
becomes controlled by T~ n \ . As shown in Fig. ||, the slope 
of the response vs. correlation plots are extremely well 
predicted by Eq. (||). 

The good agreement between the internal temperature 
selected by the system, as measured by the amplitude of 
the response of the aging system to the external per- 
turbation, and the temperature predicted theoretically 
using the IS formalism supports the main hypothesis on 
which our analysis is based, i.e. the validity of the quasi- 
equilibrium condition. The quality of the agreement is 
also a positive test on the validity of the IS formalism 
and of the proposed free energy expression. This Letter 
suggests that an aging liquid, notwithstanding its out- 
of-equilibrium condition, can still be described as a sys- 
tem in quasi-equilibrium, at the expenses of introducing 
an internal temperature which is a function both of the 
thermostat temperature and of the (^-dependent) state 
of the system — expressed by its e/s value. Eq. (0) more 
precisely defines the concept of fictive temperature, usu- 
ally defined by the experimentalist as the temperature 
at which the ergodicity of the system was broken |p6| . 
Eq. (||) clarifies that Ti nt is function also of Tf. We 
stress that the detailed analysis presented in this Let- 
ter for a structural glass is conceptually identical to the 
analysis performed in recent years to describe the out 
of equilibrium dynamics of disordered spin models p9| , 
even though we prefer to present our results in terms of 
extended validity (as opposed to violation) of the fluctu- 
ation dissipation theorem. Tests of the mean-field theory 
on finite-size disordered p-spin system also support the 



validity of the thermodynamic approach |M ■ As noted 
for the case of disordered spin-systems |2f| , the measure- 
ment of the internal temperature in aging experiments 
]3~i| may provide information, once properly interpreted, 
on dS con f I deis i.e. on structural properties of the sys- 
tem which control the equilibrium dynamics as well. Fi- 
nally, we call attention on the fact that the time window 
accessed by the numerical experiments is very different 
from the experimental one. Measurements of the internal 
temperature in aging experiments are very important to 
assess the range of validity of the presented approach. 
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FIG. 1. Basin free energy for the BMLJ system in 
harmonic approximation. Each of the dashed line shows 
fbasin (eis,T) for one specific eis value (the corresponding eis 
value coincides with the T — value of f basin)- In all curves, 
the e/s-independent contribution (3N — 3)fcsTin(fcsT) (see 
Eq.^) has not been included for clarity reasons. Filled circles 
are the equilibrium values fbasin(eis{T eq ),T eq ) . The arrows 
indicate the path followed by the system after a quench from 
Ti = 0.8 to Ti — 0.25 The inset shows the e/s-dependence of 
y^ 3 —1 3 ln(hu}i). The value of h is such that the dimension of 
Tiu) (as well as the dimension of ksT) are in units of the LJ 
potential depth m6 1. 
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10" 10 
t w (MD steps) 

FIG. 2. Left: Solutions of Eq. (|4j) for several Tf values 
for the studied BMLJ system. Right: eis as a function of 
time, following a T-jump from Ti = 0.8 to Tf = 0.25. The 
arrows show graphically the procedure which connects the 
£is(t) value to the T int value, once Tf is known. Note that 
if the curvature of the basins were independent on eis, then 
curves for different Tf would all coincide with T eq (eis) (filled 
symbols) . 



4 




FIG. 3. Left: time dependence of the responce (p^, open 
symbols) and correlation function (S^ a , filled symbols) for 
t w = 1024 (circles) and t w = 16384 (squares). Right: para- 
metric plot (in r) of p% vs S% a for the two studied t w . Dashed 
lines have slope Vo/fceT/, thick lines have slope VojksTint- 
The two Tint values can be read from FigQ The modulus of 
k is 6.7 @. 
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